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The two Higgs boson doublets in the minimal super- 
symmetric standard model (MSSM) [1] lead to five phys- 
ical Higgs bosons: three neutral (collectively denoted as 
4>): h, H, and A; and two charged: H + and H~ . Two 
parameters, conventionally chosen as the ratio of the two 
Higgs doublet vacuum expectation values, tan/3, and the 
mass of the pseudoscalar A, Ma, are sufficient to describe 
the MSSM Higgs sector at tree level. Though tan (i is a 
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free parameter in the MSSM, there are indications which 
suggest that it should be large. A value of tan j3 ~ 35 
naturally explains the top to bottom quark mass ratio [2] , 
and high tan (3 values also provide a good explanation for 
the observed density of dark matter [3]. 

The couplings of the Higgs bosons to fcrmions in the 
MSSM are proportional to the corresponding couplings 
in the standard model (SM). The proportionality factor 
depends on the type of the quark (up- or down- type) and 
on the type of the Higgs boson. At large tan (3, the two 
Higgs bosons A and either h or H have approximately the 
same mass and a down-type quark coupling enhanced by 
tan/? compared to the SM coupling, while the coupling 
to up-type quarks is suppressed. Here, the three neutral 
Higgs boson couplings to b quarks follow the sum rule 

9lbb + 9 2 Hbb + 9 Abb ~ 2 x tan2 I 3 x 9 2 h,sm, where 9H..SM 
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is the SM coupling. Therefore, in these cases the pro- 
duction of Higgs bosons associated with bottom quarks 
(down-type quarks with the highest mass) is enhanced by 
a factor of 2 x tan 2 ft compared to SM production. Due 
to the tan ft enhancement, the main decay for the three 
neutral Higgs bosons is </> — > bb with branching ratios 
near 90% (the remainder being mostly <f> — > tt) . Since a 
direct search for <f> — > bb is difficult due to large multijet 
backgrounds, searches rely on the case where <fi is pro- 
duced in association with one b quark. The final state 
with three b quarks represents a powerful search chan- 
nel, with the third 6-jet providing additional suppression 
of the large multijet background at a hadron collider. 

MSSM Higgs boson production has been studied 
at the CERN LEP e+e~ collider which excluded 
M hiA < 93 GcV for all tan/3 values [4]. The CDF [5] 
and DO [6-8] collaborations have extended MSSM Higgs 
boson searches to higher masses for high tan/3 values. 
This Letter uses data collected during Run II at the 
Fermilab Tevatron collider by the DO collaboration cor- 
responding to an integrated luminosity of 5.2 fb -1 [9]. 
The dataset is broken into two periods, corresponding 
to the period before (1.0 ftT 1 ) and after (4.2 ftT 1 ) the 
upgrade of the DO silicon vertex detector and trigger sys- 
tem. The dataset is five times larger than that used in 
the previous publication [7]. The full dataset has been 
reanalyzed to incorporate recent improvements to anal- 
ysis procedures, algorithms, and calibrations. Improved 
modeling of the background has resulted in reduced sys- 
tematic uncertainties. In addition, the Higgs boson mass 
range under consideration has been extended to 300 GeV. 
The limits are calculated using a program [10], which is 
an implementation of the modified frequentist limit set- 
ting procedure [24], and are based only on the shape, 
and not the normalisation, of the distribution of the final 
discriminating variable. 

The DO detector is described in Rcf. [11]. Dedicated 
triggers for the three trigger levels (LI, L2, L3) designed 
to select events with at least three jets are used in this 
analysis. The majority of the data were recorded with 
6-tagging requirements at the trigger level, either at L3 
or at both L2 and L3. The trigger has an efficiency of 
approximately 60% for <fib — > bbb events with a Higgs 
boson mass of 150 GeV when measured relative to events 
with three or four reconstructed jets. 

The midpoint cone algorithm [12] with a radius 1Z = 
\J (Ay) 2 + (A</?) 2 = 0.5, where y is the rapidity and ip 
the azimuthal angle, is used to reconstruct jets from 
energy deposits in the calorimeters. Details of the jet 
reconstruction and energy scale determination are de- 
scribed in Ref. [13]. In addition to passing a set of 
quality criteria, all jets are required to be matched to 
at least two tracks reconstructed in the central detec- 
tor, pointing to the pp vertex and with hits in the sili- 
con detector. The matching criterion is A7\L(track, jet- 
axis) = \J (At]) 2 + (Atp) 2 < 0.5, where 77 is the pseudo- 



rapidity. Signal events are selected by requiring three 
or four jets with transverse momenta pr > 20 GeV and 
\r)\ < 2.5. The distance in the coordinate along the beam 
axis of the position of the pp vertex (zpy) is required to 
be within 35 cm of the center of the detector. This is 
well within the geometric acceptance of the silicon detec- 
tor, as needed for efficient 6-tagging. A neural network 
(NN) 6-tagging algorithm [14], which considers lifetime 
based information involving the track impact parame- 
ters and secondary vertices, is used to identify 6-quark 
jets. Each event must have at least three jets satisfying a 
&-tag NN requirement. The single jet 6-tagging efficiency 
is approximately 50% for a light-quark jet mistag rate of 
0.8%. Data events with two tagged jets are used together 
with simulated events with two and three tagged jets to 
model the background. Finally, the transverse momenta 
of the two 6-tagged jets with the highest pr are required 
to be above 25 GeV. The data are split into four inde- 
pendent channels based on jet multiplicity (three or four 
jets) and running period. 

The leading order event generator pythia [15] is used 
to generate samples of associated production of <j> and a 
b quark in the 5-flavor scheme [16], gb — > <f)b. The cross 
section [17], experimental acceptance, and the kinematic 
distributions of the 6-quark jet produced in association 
with the Higgs boson are corrected to next-to-leading or- 
der (NLO) using mcfm [16]. Multijet background events 
from the bbj, bbjj, ccj, ccjj, bbec, and bbbb processes, 
where j denotes a light parton (u, d, s quark or gluon), 
are generated with the ALPGEN [18] event generator. A 
matching algorithm is used to avoid double counting of fi- 
nal states [19]. The small contribution from production 
to the background is also simulated with ALPGEN. Other 
processes, such as Zbb and single top quark produc- 
tion, are negligible. The ALPGEN samples are processed 
through pythia for showering and hadronization. All 
samples are then processed through a GEANT-based [20] 
simulation of the DO detector. The same reconstruction 
algorithms are used for the simulated samples as for the 
data. A parameterized trigger simulation, based on ef- 
ficiencies measured in data, is used to model the effects 
of the trigger requirements. The 6-tagging is modeled by 
weighting simulated events based on their tagging prob- 
ability measured using data [14]. The efficiency of the 
requirements on the triggers, zpy, and the number of 
jets, range from 1.9% to 26.4% for Higgs boson masses 
between 90 and 300 GeV. After the three 6-tag require- 
ment the efficiency with respect to the total number of 
signal events ranges from 0.2% to 1.4% in the three-jet 
channel (0.1% to 0.9% in the four-jet channel). Table I 
shows the number of events in data and the signal effi- 
ciency at different stages of the event selection. 

Multijet processes contribute to the background and 
the theoretical uncertainty on their cross sections is very 
large. The background composition is therefore deter- 
mined by fitting distributions of the transverse momenta 
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Number of events Fraction (%) 


Signal eff. (%) 




(xlO 3 ) 




M A = 200 GeV 


Events 


517,288 


100 


- 


Trigger 


198,106 


38 


30 


Zpv cut 


195,587 


38 


25 


3/4 jets 


96,318/21,898 


19/4.2 


12/3.7 


2 6-tag jets 


710/230 


0.14/0.044 


5.4/1.8 


3 6-tag jets 


15/11 


0.0029/0.0021 


1.2/0.61 



TABLE I: The number and fraction of events in data and 
signal efficiency for each selection requirement. As the data 
are split into three- and four-jet sub-samples, these numbers 
are reported separately in the last three rows. 



of jets of simulated events to data. The fractional contri- 
bution, on, of the ith multijet background process is cal- 
culated from equations linking the 6-tag efficiency for the 
ith background, e\, with the number of observed events, 
Nk, where k indicates the number of 6-tagged jets (0-3) 
in an event [21], and the total number of events, N tot : 



J2i x e l 



N k /N u 



(1) 



The double 6-tagged sample is dominated by bbj, while 
the triple 6-tagged sample consists of a mixture of approx- 
imately 50% 666, 30% bbj, 15% bbc+bec and a remaining 
fraction consisting of ccj, bjj, cjj, and jjj. 

For every event, the two jet pairs with the largest scalar 
summed transverse momenta are considered as possible 
Higgs boson candidates. To remove discrepancies be- 
tween data and simulation originating from gluon split- 
ting (g — > 66), jet pairs which do not fulfill A1Z > 1.0 are 
rejected. 

Six variables for which the data distributions are well 
modeled by the simulation are used to separate the jet 
pair from a Higgs boson from the background: At] be- 
tween the two jets in the pair, A<j> between the two jets in 
the pair, the angle between the leading jet in the pair and 
the total momentum of the pair, the momentum balance 
in the pair [22], the combined rapidity of the pair, and 
the event sphericity. Based on these kinematic variables 
a likelihood discriminant, V, is calculated according to: 



V(x 1 ,....,x e ) 



nli^ SIS (^) 



n-=i^ Ig (^)+nr =1 ^ bkg (^) 



(2) 



where P^ 1S (Pi ks ) refers to the signal (background) 
probability density function (pdf) for variable Xi, and 
(xi, xg) is the set of measured kinematic variables. 
The pdfs are obtained from triple 6-tagged signal and 
simulated background samples. Two likelihood discrimi- 
nants, providing discrimination in the Higgs boson mass 
ranges 90 - 130 GeV (low-mass) and 130 - 300 GeV 



(high-mass), respectively, are built by combining simu- 
lated signal samples from the appropriate Higgs boson 
mass ranges. Signal samples of equal size, interspaced 
by 10 GeV in Ma, are hence added together within the 
low-mass and high-mass range, respectively. After eval- 
uating the likelihood, only the jet pairing with the larger 
T> is kept for each event in each mass range. To fur- 
ther remove background from the final analysis sample, 
events are only selected if V > 0.65. The likelihood re- 
quirements are optimized considering the variation of the 
predicted limit in tan /3. The final discriminant used in 
the limit calculation is the distribution of the jet pair in- 
variant mass, M bb , after the selection requirement of the 
likelihood appropriate to the mass of the hypothesized 
Higgs boson. 

The 666 background is indistinguishable from the sig- 
nal events where the wrong 6-jet pair is chosen by the 
likelihood and consequently cannot be normalised from 
the data. The M bb background shape is modeled using 
a combination of data and simulated samples. The dis- 
tribution, Sy^^ g {V, M bb ), of the predicted triple 6-tagged 
(3Tag) background sample in the two-dimensional T> and 
M hb plane is obtained from the inclusive double 6-tagged 
(2Tag) data shape multiplied by the ratio of the simu- 
lated shapes of the SM triple (5*|xag) an d double tagged 
events (S™%): 



oMC Ijy tut _\ 
oPrcd (T) M _\ _ 3TagV > bbj „ data „ _x 

S 3Tae (V,M bb ) - s ^ c Mb _ b) i> 2Tag (V,M bb ). 



(3) 



Fig. 1 shows T> for data and background for the low- and 
high-mass likelihoods in the three-jet channel. The shape 
of a signal distribution, normalised to the same number 
of events as data, is also shown. Fig. 2 shows the M bb 
distribution in the three-jet channel after the low- and 
high-mass likelihood selection requirements, respectively. 
The invariant mass of a Higgs boson signal in the three- 
jet channel is shown in Fig. 3 for three different values of 
M A . 

To verify the background model a signal-depleted re- 
gion is studied - any deviation observed there is unlikely 
to be as a result of signal and therefore would indicate a 
possible problem in the background modeling. A sample 
is hence chosen using the lower likelihood jet pairing and 
applying a selection of V < 0.12. Fig. 4 shows the invari- 
ant mass distributions for background and data in this 
sample. Agreement (x 2 /n.d.f. = 0.86) between the back- 
ground model and the data is observed. A wide variety 
of additional cross-checks were carried out, examining as- 
pects of the event selection, b-tagging, and background 
modelling; no significant changes in the results were ob- 
served. 

Sources of systematic uncertainty on both the signal 
normalisation and shape are considered. The sources 
of systematic uncertainty on the signal included are: 6- 
quark jet identification efficiency, 6- and light-quark jet 
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FIG. 1: (color online) Comparison of the a) low-mass and b) high-mass likelihood distributions for the data and predicted 
background defined by Eq. 3 in the 3Tag exclusive three-jet channel. Each event has one entry, the jet pairing with the highest 
likelihood output. Black crosses refer to data, the solid line shows the total background estimate, and the shaded region 
represents the heavy flavor component (bbb, bbc, and bcc). The distributions for a Higgs boson of mass 120 and 200 GeV are 
shown as a dashed line in a) and b), respectively. The background and signal contributions are normalised to have an equal 
number of events as data. The arrows indicate the selection cut at V — 0.65. 



energy resolution, trigger modeling, jet energy calibra- 
tion, jet identification, integrated luminosity, and theo- 
retical models. The theoretical uncertainty on the signal 
cross section is estimated from MCFM [16] and consists 
of a contribution of 10% from the choice of factorisation 
scale as well as an uncertainty of 5% to 13% from the 
parton distribution functions, depending on Higgs boson 
mass. Both the theoretical uncertainty and the luminos- 
ity uncertainty of 6.1% [9] are treated as normalisation 
uncertainties for each mass hypothesis. The remaining 
sources of systematic uncertainty are assessed individu- 
ally by varying parameters within their uncertainties and 
taking into account the resulting difference in normali- 
sation and shape of the M bb distribution at each mass 
point. 

For the dominant background, only systematic uncer- 
tainties affecting the shape of M bb matter, since only 
the shape and not the normalisation is used to distin- 
guish signal from background in this analysis. Addition- 
ally, many uncertainties affecting the simulation, like the 
jet energy scale and resolution uncertainties, cancel in 
Eq.3. The estimated variations in the remaining system- 
atic sources are propagated to T>, M bb and the predicted 
shape S'fxlg an d used in assessing the limits presented 
below. The uncertainty from the ^-tagging of jets is eval- 
uated by varying the 6-tag efficiencies within their uncer- 
tainties. The uncertainties in the difference of the energy 
resolution between heavy and light flavor jets is obtained 
by smearing the energy of the b- and c-quark jets by an 
additional 7%, corresponding to half the light-quark jet 



energy resolution. The shape difference between triple 
and double 6-tagged data in the trigger turn-on curves 
resulting from the 6-tagging criteria in the trigger is ac- 
counted for as a systematic uncertainty. Small variations 
in the shape, arising from possible signal contamination 
when determining the background composition, are in- 
cluded as a systematic uncertainty. Finally, the uncer- 
tainty on the ti normalisation is taken as 10% [23]. 

No significant excess over the background is observed 
in the data. Limits on the Higgs boson production 
cross section multiplied by the branching ratio to bb 
are therefore calculated with the modified frequentist 
method [10, 24]. The confidence level of the signal, CL S , 
which is used to calculate the exclusion, is defined in 
Ref. [24]. The overall normalisation of the background 
expectation is allowed to float independently in the null 
(background-only) and test (background-plus-signal) hy- 
potheses. The systematic uncertainties on the signal and 
background are included in the limit setting procedure. 
Each component of systematic uncertainty is adjusted 
by introducing multiplicative scale factors and maximiz- 
ing the likelihood for the agreement between prediction 
and data with respect to these scale factors, constrained 
by prior Gaussian uncertainties. Limits on the prod- 
uct of cross section and branching ratio are obtained by 
scaling the signal cross section until 1 — CL S = 0.95 is 
reached. These limits are effectively independent of the 
signal model but assume the width of the <f> to be neg- 
ligible relative to the experimental resolution («20% at 
Ma = 150 GeV). The four independent analysis chan- 



7 



> 

<D 

(5 
o 

■*— < 
c 

> 

LU 



D0 , 5.2 fb 1 
a) 3 jet 

Low-mass likelihood 

—I— D0 Data 

Background 

Heavy flavor 




^1200 

|iooo 

I 800 
ui 



600 

400 

200 


100 



D0 , 5.2 ft> 1 
b) 3 jet 

High-mass likelihood 

— I— D0 Data 

Background 

Heavy flavor 




50 100 150 200 250 300 350 400 

M bB [GeV] 



50 100 150 200 250 300 350 400 

M bB [GeV] 



FIG. 2: (color online) Invariant mass distribution for the 3Tag exclusive three-jet channel for: a) the low-mass likelihood 
selection and b) the high-mass likelihood selection. Each event has one entry, the jet pairing with the highest likelihood output. 
Black crosses refer to data, the solid line shows the total background estimate, and the shaded region represents the heavy 
flavor component (bbb, bbc, and bcc). The lower panels show the difference between the data and the predicted background. 



nels are combined in the limit setting procedure. Signal 
hypotheses are considered for discrete Higgs boson mass 
points from 90 to 300 GeV in steps of 10 GeV. The treat- 
ment of the systematic uncertainties and the limit setting 
procedure were extensively cross-checked; no unexpected 
effects were observed. 

The combined result is summarized numerically in Ta- 
ble II and the model independent limit is shown in Fig. 5. 
The deviation from expectation around 120 GeV corre- 
sponds to 2.5 standard deviations. Note that it is more 
likely to find a deviation (in the background-only hypoth- 
esis) when several mass bins are probed than if only one 
bin is probed. A standard convention to account for this 
"trial factor" [25] gives a significance of the deviation at 
120 GeV of 2.0 standard deviations. 

As a consequence of the enhanced couplings to 6-quarks 
at large tan/3, the total width of the Higgs boson mass 
also increases with tan f3. This can have an impact on the 
search if the width is comparable to or larger than the 
experimental resolution. To take this effect into account, 
the width of the Higgs boson is calculated with feyn- 
HIGGS [26] and included in the simulation as a function 
of the mass and tan /3 by convoluting a relativistic Breit- 



Wigner function with the NLO cross section. The masses 
and couplings of the Higgs bosons in the MSSM depend, 
in addition to tan/3 and Ma, on the SUSY parameters 
through radiative corrections. Limits on tan/3 as a func- 
tion of Ma are derived for two particular scenarios as- 
suming a CP-conserving Higgs sector [27]: the m™ ax [28] 
and no-mixing [29] scenarios with a negative or positive 
value of the Higgs sector bilinear coupling, fi. Figure 6 
shows the result interpreted for these two scenarios in the 
case of /i = —200 GeV. Weaker limits are obtained for 
the fi > scenarios, due to the decrease in the product 
of cross section and branching ratio for positive values of 
/i [27]. 

The results exclude substantial areas in the MSSM pa- 
rameter space up to Higgs boson masses of 300 GeV, 
under the assumption that a perturbative treatment is 
valid over the entire region. These are the most strin- 
gent limits to date for this topology over this mass range 
at a hadron collider. 

We thank the staffs at Fermilab and collaborating 
institutions, and acknowledge support from the DOE 
and NSF (USA); CEA and CNRS/IN2P3 (France); 
FASI, Rosatom and RFBR (Russia); CNPq, FAPERJ, 
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signals in the three-jet channel for Ma = 120, 180 and 240 
GeV, respectively. The distributions are normalised to unit 
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TABLE II: Model independent 95% C.L. upper limits on the 
cross section times branching ratio for the combined 5.2 fb _1 
analysis. 
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FIG. 4: (color online) Invariant mass distribution for the jet 
pair with the lower likelihood for data (black crosses) and 
background (line) in the 3Tag exclusive three-jet channel for 
the 4.2 fb _1 data set collected after the silicon detector up- 
grade. The lower panel shows the difference between the data 
and the predicted background. 



search Council (Sweden); and CAS and CNSF (China). 



[1] 

[2] 

[3] 
[4] 

[5] 

[6] 
[7] 
[8] 



H.P. Nilles, Phys. Rep. 110, 1 (1984); H.E. Haber and 
G.L. Kane, Phys. Rep. 117, 75 (1985). 
B. Ananthanarayan, G. Lazarides, and Q. Shafi, Phys. 
Rev. D 44, 1613 (1991). 

V. Barger and C. Kao, Phys. Lett. B 518, 117 (2001). 
The ALEPH Collaboration, The DELPHI Collaboration, 
The L3 Collaboration, and The OPAL Collaboration, 
Eur. Phys. J. C 47, 547 (2006). 

T. Affolder et al. (CDF Collaboration), Phys. Rev. Lett. 
86, 4472 (2001); A. Abulencia et al. (CDF Collabora- 
tion), Phys. Rev. Lett. 96, 011802 (2006). 
V.M. Abazov et al. (DO Collaboration), Phys. Rev. Lett. 
95, 151801 (2005). 

V.M. Abazov et al. (DO Collaboration), Phys. Rev. Lett. 
101, 221802 (2008). 

V.M. Abazov et al. (DO Collaboration), Phys. Rev. Lett. 
101, 071804 (2008); V.M. Abazov et al. (DO Collabora- 
tion), Phys. Rev. Lett. 102, 051804 (2009); V.M. Abazov 
et al. (DO Collaboration), Phys. Rev. Lett. 104, 151801 
(2010). 



9 



D0, 5.2fb 



Expected 
Expected ±1 s.d. 
Expected ±2 s.d. 
— Observed 




100 



150 



200 



250 300 
M A [GeV] 



FIG. 5: (color online) Model independent 95% C.L. upper 
limit on the cross section multiplied by the branching ratio 
for the combined 5.2 fb _1 analysis. The light and dark grey 
regions correspond to the one and two standard deviation 
(s.d.) bands around the expected limit. 



[9] T. Andeen et al., FERMILAB-TM-2365 (2007). 
[10] W. Fisher, FERMILAB-TM-2386-E (2007). 
[11] V.M. Abazov et al. (DO Collaboration), Nucl. Instrum. 

Methods Phys. Res. A 565, 463 (2006); V.M. Abazov et 

al. (DO Collaboration), Nucl. Instrum. Methods in Phys. 

Res. Sect. A 622, 298 (2010). 
[12] G. Blazey et al, arXiv:hep-ex/0005012 (2000). 
[13] V.M. Abazov et al. (DO Collaboration), Phys. Rev. Lett. 

101, 062001 (2008). 
[14] V.M. Abazov et al. (DO Collaboration), Nucl. Instrum. 

Methods Phys. Res. A 620, 490 (2010). 
[15] T. Sjostrand, S. Mrenna, and P. Skands, J. High Energy 



Phys. 05, 026 (2006). Version 6.409 is used. 

[16] J. Campbell, R.K. Ellis, F. Maltoni, and S. Willenbrock, 
Phys. Rev. D 67, 095002 (2003). 

[17] All cross sections are quoted with respect to acceptance 
requirements on the highest-pT 6-quark not from the 
Higgs decay with pr > 12 GeV and \r]\ < 5.0. 

[18] M.L. Mangano et al, J. High Energy Phys. 07, 001 
(2003). ALPGEN version 2.11 is used. 

[19] S. Hoche et al, arXiv:hep-ph/0602031 (2006). 

[20] R. Brun and F. Carminati, CERN program library long 
writeup W5013 (1993). 

[21] The events with four 6-tagged jets in the four jet sam- 
ple are included with the events with three 6-tagged jets 
when determining the background composition. 

[22] The momentum balance is defined as \ph 1 — Pb 2 \/\Pb 1 + 
pt 2 1, where is the magnitude of the momentum three 
vector of the ith 6-quark jet. 

[23] S. Moch and P. Uwer, Nucl. Phys. Proc. Suppl. 183, 75 
(2008); S. Moch and P. Uwer, Phys. Rev. D 78, 034003 
(2008). 

[24] T. Junk, Nucl. Instrum. Methods Phys. Res. A 434, 435 

(1999) ; A. Read, Nucl. Instrum. Methods Phys. Res. A 
425, 357 (1999). 

[25] The ALEPH Collaboration, The DELPHI Collaboration, 
The L3 Collaboration, The OPAL Collaboration, and 
The LEP Working Group for Higgs Boson Searches, Phys 
Lett B 565, 61 (2003). 

[26] S. Heinemeyer, W. Hollik, and G. Weiglein, Eur. Phys. 
J. C 9, 343 (1999); Comput. Phys. Commun. 124, 76 

(2000) ; G. Degrassi et al, Eur. Phys. J. C 28, 133 (2003); 
M. Frank et al, J. High Energy Phys. 02, 047 (2007). 
feynhiggs version 2.6.5 is used. 

[27] M. Carena, S. Heinemeyer, C. E. M. Wagner, and 

G. Weiglein, Eur. Phys. J. C 45, 797 (2006). 
[28] Msusy = 1 TeV, X t = 2 TeV, M 2 = 0.2 TeV |/u| = 0.2 



TeV, and m„ 



0.8 TeV. 



[29] M s 



USY 



2 TeV, X t = TeV, M 2 = 0.2 TeV \fi\ = 0.2 



TeV, and m„ 



1.6 TeV. 



10 




FIG. 6: (color online) a) 95% C.L. lower limit in the (M^,tan/3) plane obtained for the m™ ax , /i = —200 GeV scenario, b) the 
lower limit for the no-mixing, n — —200 GeV scenario. The one and two standard deviation bands around the expected limit 
and the exclusion limit obtained from the LEP experiments are also shown [4]. 



